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from  near-field  measurements  of  upholo  time  and  spall  observations  for  47  ex- 
plosions. Spall  energies  relative  to  explosion  energies  were  estimated  as- 
suming that  a single  spall  closure  was  responsible  for  most  observed  surface 
displacements.  Some  calculated  spall  energies  were  larger  than  the  explosion 
energy  entering  the  spall  region,  suggesting  that  the  method  for  calculating 
spall  energy  is  in  error.  Relative  spall  energies  compared  to  explosion 
energies  foi  42  NTS  explosions  located  below  the  water  table  show  no  relation- 
ship to  Mg  residuals  or  strain  release.  Subsurface  vertical  accelerometer  re- 
cords show  that  spall  energy  is  divided  among  several  closures  over  a range  of 
depths,  the  deeper  formations  closing  before  the  shallower  layers.  The  large 
observed  spall  pulses  are  due  to  closures  closest  to  the  gauge.  The  principal 
error  in  the  spall  energy  calculat ions  was,  therefore,  the  assumption  that  the 
first  spall  closure  contained  most  of  the  spall  energy.  The  total  spall 
energy,  the  sum  of  many  spall  closures,  is  probably  much  smaller  than  the  spall 
energy  calculated  for  a single  spall  closure.  The  subsurface  accelerometer 
data  show  little,  if  any,  spall  energy  leaving  the  source  region,  suggesting 
that  spall  has  little  effect  on  and  Ms  at  teleseismic  distances. 
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This  study  examines  near-field  observations  of  spalling,  a phenomenon 
associated  with  most  underground  nuclear  explosions,  and  attempts  to  esti- 
mate its  effects  on  m^  and  Ms.  Observations  of  spall  on  near-field  accel- 
erometers and  particle  velocity  records  were  used  to  examine  possible  re- 
lationships between  radius  of  spallation,  estimated  spall  thickness,  max- 
imum height  attained  by  spall  at  the  surface,  shot  depth  of  burial,  and 
yield. 

Estimates  of  pP-P  and  Pe_.rr-P  at  teleseismic  distances  were  calculated 
from  near-field  measurements  of  uphole  time  and  spall  observations  for  47 
explosions.  Spall  energies  relative  to  explosion  energies  were  estimated 
assuming  that  a single  spall  closure  was  responsible  for  most  of  the  observ- 
ed surface  displacements.  Some  calculated  spall  energies  were  larger  than 
the  explosion  energy  entering  the  spall  region,  suggesting  that  the  method 
for  calculating  spall  energy  is  in  error.  Relative  spall  energies  compared 
to  explosion  energies  for  42  NTS  explosions  located  below  the  water  table 
show  no  relationship  to  residuals  or  strain  release.  Subsurface  verti- 
cal accelerometer  records  show  that  spall  energy  is  divided  among  several 
closures  over  a range  of  depths,  the  deeper  formations  closing  before  the 
shallower  layers.  The  large  observed  spall  pulses  are  due  to  closures 
closest  to  the  gauge.  The  principal  error  in  the  spall  energy  calculations 
was,  therefore,  the  assumption  that  the  first  spall  closure  contained  most 
of  the  spall  energy.  The  total  spall  energy,  the  sum  of  many  spall  closures, 
is  probably  much  smaller  than  tht  ;pall  energy  calculated  for  a single  spall 
closure.  The  subsurface  accelerometer  data  show  little,  if  any,  spall 
energy  leaving  the  source  region,  suggesting  that  spAll  has  little  effect 
on  m^  and  at  teleseismic  distances. 


TABLE  OF  CONTENTS 


ABSTRACT 
LIST  OF  FIGURES 
LIST  OF  TABLES 
INTRODUCTION 

NEAR-FIELD  MEASUREMENTS  OF  SPALL 
ESTIMATES  OF  SPALL  PARAMETERS 
ESTIMATES  OF  TELESEISMIC  pP  AND  P 

o PALjL# 

ARRIVAL  TIMES  FROM  NEAR-FIELD  DATA 
SPALL  ENERGY  AND  ITS  EFFECTS  ON  AND  Mg 
CONCLUSIONS 
ACKNOWLEDGEMENTS 
REFERENCES 

APPENDIX  - Surface  Particle  Velocity  Data 


Page 

2 

4 

5 

7 

8 

16 

32 

35 

51 

54 

A-l 


LIST  OF  FIGURES 


F igure 
1 
2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


No.  Title 

Generation  of  a spall  by  a compressional  stress  wave. 

Direct  ray  and  spall  pulse  recorded  at  the  surface  near 
an  explosion. 

Variation  of  maximum  height  of  projection  of  spall  layer  (s) 
vs.  range  (5). 

Surface  accelerometer,  particle  velocity  and  displacement 
data  at  various  horizontal  distances  from  the  shot  hole 
for  the  Aleutian  Islands  explosion  M1LR0W. 

Thickness  of  spall  vs.  yield. 

Effective  spall  thickness  (corrected  to  122  mkt  ) vs. 
yield. 

Scaled  spall  thickness  vs.  scaled  shot  depth  of  burial. 

Spall  thickness  vs.  shot  depth  of  burial. 

Radius  of  spallation  vs.  yield. 

Scaled  spall  radius  vs.  scaled  shot  depth  of  burial. 

S (maximum  height  attained  by  spall)  observed  near  the 
source  vs.  scaled  shot  depth  of  burial. 

2 

S calculated  from  (At)  g/8  vs.  scaled  shot  depth  of 

.max  , 

burial. 

S observed  near  the  source  vs.  yield, 
max 

Spall  radius  vs.  spall  thickness. 

Energy  vs.  range  for  five  spall  closures. 

Yield  vs.  radius  of  maximum  spall  energy. 

Energy  vs.  time  for  spall  closure. 

Spall  energy /explosion  energy  vs.  yield. 

M residual  vs.  spall  energy/explosion  energy  for  A2  NTS 
explosions  located  below  the  water  table. 


19 


LIST  OF  FIGURES  (continued) 


Figure  No.  Title  Page 

20  Subsurface  vertical  accelerometer  array  data  for  RAINIER  46 

(Eisler,  Chilton  and  Sauer,  1966). 

21  Subsurface  vertical  accelerometer  array  data  for  GNOME  47 

(Weart,  1962). 

22  Subsurface  vertical  accelerometer  array  data  for  SALMON  48 

(Eisler , 1967) . 

23  Subsurface  vertical  accelerometer  array  data  for  MILROW  49 

(Perret  and  Breding,  1972) . 


- 5 - 


LIST  OF  TABLES 


Table  No. 

Title 

Page 

I 

Spall  parameters. 

29 

II 

Teleseismic  pP-P  and  P„  -P  estimates. 

SPALL 

33 

III 

Spall  energies. 

40 

IV 

Explosion  energies  from  Perrot  (1972)  and  from  Perret's 
efficiency  estimates  for  the  explosion  medium. 

44 

r 


INTRODUCTION 

Spalling  is  a phenomenon  associated  with  most  underground  nuclear  explo- 
sions. The  compressional  stress  wave  leaving  the  source  is  reflected  at  the 
Earth's  surface  and  travels  downward  as  a tensile  stress  wave.  Where  tensile 
stress  exceeds  compressive  stress  remaining  in  the  upward  travelling  wave, 
plus  stress  due  to  the  overburden  and  tensile  strength  of  the  rock,  near- 
surface layers  will  physically  part  (Figure  1) . The  impact  during  closure  of 
the  spall  gap  produces  a signal  observed  on  near-field  instruments.  The 
purpose  of  this  study  is  to  estimate  the  effect,  if  any,  of  spall  on  tele- 
seismic  and  measurements. 

Estimated  spall  depth,  spall  radius,  maximum  height  attained  by  spall, 
shot  depth  of  burial,  and  yield  were  examined  to  determine  any  possible  re- 
lationships between  these  parameters.  Estimates  of  teleseismic  Pcr,ATT_p 

brALl.* 

times  were  calculated  from  near-field  accelerometer  data  for  47  explosions 
with  observed  spall.  Spall  and  explosion  energies  were  calculated  and  sub- 
surface vertical  accelerometer  data  examined  to  determine  the  amount  of  spall 
energy  leaving  the  source  region.  This  study  demonstrates  that  assuming  the 
first,  and  deepest,  spall  closure  contains  most  of  the  spall  energy  is  wrong 
and  that  calculations  of  spall  energy  cannot  be  based  solely  upon  near-field 
surface  measurements. 


NEAR-FIELD  MEASUREMENTS  OF  SPALL 


't 


An  underground  explosion  produces  spalling  when  the  tensile  stress  wave 
reflected  at  the  surface  exceeds  compressive  stress  remaining  in  the  upward 
travelling  wave,  plus  the  stress  due  to  the  overburden  and  the  tensile  strength 
of  the  rock  (Figure  1).  The  energy  of  the  waves  trapped  in  the  formation 
causes  the  layers  to  ascent]  and  spall  is  created.  At  the  surface  near  an 
explosion,  the  first  signal  recorded  is  the  pulse  travelling  upward  from  the 
explosion,  t^.  The  second  signal,  t^,  originates  from  the  spall  closure 
directly  below  the  point  of  observation.  Figure  2 shows  the  two  ray  paths 
schematically.  The  arrivel  time  between  t and  t is  denoted  At.  The  first 
signal  on  the  accelerometer  trace  is  a positive  pulse  from  the  explosion  fol- 
lowed by  a free  fall  period  of  -g,  and  then  a positive  pulse  from  the  spall 
closure. 

In  several  already  existing  studies,  investigators  assume  a constant  spall 
layer  thickness  since  no  observations  of  lateral  variations  in  spall  thickness 
exist.  The  maximum  height  the  layer  is  projected,  is  generally  observed 

to  be  greater  over  the  explosion,  tapering  off  to  zero  at  increasing  distances 
from  the  source.  The  shape  of  spall  height  is  assumed  to  be  conical  with  sides 
concave  upward  as  shown  in  Figure  3 (Eisler  and  Chilton,  1964). 

Since  the  maximum  observed  surface  displacement  is  greater  than  the  S 

max 

derived  from  displacement  data  both  above  and  below  the  spall  depth,  addition- 
al spalling  at  other  depths  must  bo  taking  place  (Eisler  et.  al.,  1966). 

Several  possible  spall  signals  are  recorded  at  the  surface  for  some  events 
and  subsurface  near-field  instrumentation  indicates  several  spall  depths  for 
some  events.  After  the  first  spall  is  created,  remaining  parts  of  the  stress 
wave  may  travel  upward  and  downward  to  create  more  spalls  above  and  below  the 
first  spall.  The  subsurface  spalled  region  may  have  upper  and  lower  bounds 
where  the  pulses  erode  to  the  point  that  spalling  no  longer  takes 

Eisler,  J. , and  Chilton,  F.,  1964  Spalling  of  the  Earth's  surface  by  under- 
ground nuclear  explosions,  J.  Geophys  Res . , 69,  5285. 

Eisler,  J,. , Chilton,  F.,  and  F.  Sauer,  1966.  Multiple  subsurface  spalling 
by  underground  nuclear  explosions,  J.  Geophys . Res.,  71,  3923. 
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Figure  2.  Direct  ray  and  spall  pulse  recorded  at 
the  surface  near  an  explosion. 
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place  (Chilton  et.  al.,  1966).  The  thicknesses  of  all  the  layers  are  assumed 
constant  because  no  observations  of  their  shapes  exist. 

Figure  4 shows  several  surface  records  for  the  explosion  MILROW  in  the 
Aleutian  Islands.  Spalling  signals  were  recorded  out  to  10,000  meters  of 
horizontal  range  from  the  source.  Vertical  arrays  of  accelerometers  showed 
zones  of  spall  partings  at  100  to  300  feet,  300  to  500  feet,  and  500  to  1000 
feet.  A computer  program  that  simulates  velocities  and  stress  profiles  with 
time  in  the  rock  column  indicated  that  a parting  would  occur  at  from  350  to 
400  feet  (Perret  and  Breding,  1972)  The  spall  at  this  depth  range  was  the 
first  to  close  and  produced  the  first  and  largest  spall  signal  on  the  accel- 
erometers. Note  that  spall  signals  on  accelerometer  records  are  relatively 
large  compared  to  the  main  phase  from  the  explosion.  The  large  signals  occur 
because  the  source  of  the  spall  signal  is  directly  under  the  accelerometer  and 
it  has  travelled  a smaller  distance  and  is  less  attenuated  t) an  the  explosion 
signal  . 


Perret,  W.,  and  D.  Breding,  1972.  Ground  motion  in  the  vicinity  of  an  under- 
ground nuclear  explosion  in  the  Aleutian  Islands:  MILROW  event,  Sandia 
Laboratories,  Research  Report  SC- RR- 7 1-0668. 
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Figure  4.  (cont.)  Surface  accelerometer,  particle  velocity. 


and  displacement  data  at  various  horizontal 
distances  from  the  shot  hole  for  the  Aleutian 
Islands  explosion  MILROW. 
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Figure  4.  (cont.)  Surface  accelerometer,  particle  velocity, 

and  displacement  data  at  various  horizontal 
distances  from  the  shot  hole  for  the  Aleutian 
Islands  explosion  MILROW. 
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Viecelli  (1973)  estimated  the  depth,  height,  and  radial  extent  of 
spallation  from  near-field  accelerometer  records  for  seven  NTS  explosions 
and  plotted  these  parameters  against  shot  depth  of  burial  and  yield.  The 
present  study  will  add  more  data  points  to  his  figures  from  47  explosions 
at  NTS , as  well  as  other  test  areas  in  Nevada,  New  Mexico,  Mississippi  and 
the  Aleutians.  These  estimates  of  spall  parameters  will  be  used  in  a later 
section  to  determine  the  energy  of  the  spall  closure  for  a falling  disc  with 
parameters  estimated  from  near-field  surface  accelerometer  and  particle 
velocity  data. 

The  thickness  of  spall  was  determined  from  subsurface  vertical  acceler- 
ometer array  data  near  the  sources,  observations  of  partings  in  boreholes, 
and  determination  of  the  arrival  of  the  phase  t^,  which  is  reflected  off  the 
spall  layer  (Figure  2).  The  first  spall  closure  was  assumed  to  contain  most 
of  the  energy.  The  maximum  height  of  the  spall  was  observed  by  near-field 
vertical  arrays  or  calculated  theoretically  from  the  observed  maximum  value 
of  At  = t -tg  at  the  surface  and  the  acceleration  of  gravity  g according  to 
vhe  formula 


S 

max 


(At)  "g 
8 


(Kisler  and  Chilton,  1964).  The  radial  extent  of  spallation  is  defined  as 
that  horizontal  range  from  the  source  where  At  becomes  constant.  For  some 
explosions  this  condition  was  approached  near  the  largest  range  measured  and 
the  radius  of  spallation  is  taken  as  the  maximum  horizontal  range  of  measure- 
ment. 


i 


4 


Viecelli,  J.  , 1973.  Spallation  and  the  generation  of  surface  waves  by  an 
underground  explosion,  J.  Geophys.  Res.,  78 , 2475. 
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In  a few  cases  near-field  vertical  accelerometer  data  was  available 
from  Eisler  (1967),  Eisler  and  Chilton  (1964),  Eisler,  Chilton  and  Sauer 


(1966) , Perret  (1968) , Perret  and  Breding  (1972) , Viecelli  (1973) , and 
Weart  (1962). 


•t 


Eisler,  J.,  1967.  Near-field  spalling  from  a nuclear  explosion  in  a salt 
dome,  J.  Geophys.  Res.,  72,  1751. 

Perret,  W.,  1968.  Free-filed  particle  motion  from  a nuclear  explosion  in 

salt.  Part  I,  Vela  Uniform  Program,  SALMON  event,  Sandia  Laboratories, 
Report  VUF-3012. 

Weart,  W. , 1962.  Particle  motion  near  a nuclear  detonation  in  halite. 


In  all  cases  Sandla  Laboratory  near-field  particle  velocity  data,  that  D. 
Springer  made  available  to  the  author  at  the  Lawrence  Livermore  Laboratory, 
was  used  to  determine  values  of  At  and  ground  displacement  as  a function  of 
horizontal  range.  The  Appendix  includes  this  data  for  all  horizontal  ranges 
out  to  the  radial  extent  of  spallation.  The  surface  particle  velocity  data 
listed  times  between  the  first  and  second  zero  crossing  on  the  particle 
velocity  records  at  several  horizontal  ranges  for  each  explosion.  When  the 
particle  velocity  records  were  compared  to  the  accelerometer  data  (Figure  4) , 
these  times  were  found  to  be  equal  to  0.5At.  Note  in  the  appended  material 
that  peak  displacements  do  not  always  decrease  with  increasing  range  for  any 
given  explosion.  The  conical  approximation  for  the  shape  of  is  not 

always  valid,  probably  because  of  changing  lithologies  in  the  source  region. 

In  most  cases  this  study  picked  a larger  spall  radius  for  the  explosions 
thati  the  radius  used  in  Viecelli's  study.  This  discrepancy  stems  from  the 
methods  used  in  the  two  different  studies  to  determine  spall  radius. 

Viecelli,  as  shown  in  Figure  2,  uses  the  horizontal  radius  of  the  minimum 
arrival  time  of  t2  as  the  spall  radius.  He  reasoned  that  where  the  t2  arrival 
time  was  a minimum,  the  free  fall  time  of  the  slab  under  that  point  was 
smallest.  The  Sandia  Laboratory  particle  velocity  data  gave  values  of  At 
and  not  values  of  t2  and  tQ  at  various  ranges  from  the  source.  The  spall 
radius  was  established  as  the  range  where  the  time  At  became  constant, 
assuming  that  beyond  that  range  there  was  no  spall  under  those  observation 
points . 

Figure  5 shows  the  thickness  of  spall  versus  explosion  yield.  In  Figures 
5 through  14  the  dots  are  Viecelli's  data  points  and  the  circles  are  from 
this  study.  Data  points  for  this  study  are  tabulated  in  Table  I.  Arrows 
connect  the  two  data  points  where  our  data  points  differ  from  Viecelli's  for 
the  same  explosions.  The  figures  show  the  lines  Viecelli  used  as  best  fits 
to  his  data  and  the  least  squares  lines  of  best  fit  calculated  when  the  new 
data  points  are  added  to  the  figures.  The  additional  data  points  show  more 
scatter  than  Viecelli's  original  set,  a conclusion  found  in  all  the  figures 
resulting  from  changing  lithologies  in  the  source  regions.  Viecelli's  line 
in  Figure  5 has  a slope  of  1/3.  The  present  study  also  indicates  the  same 
relationship  between  spall  thickness  and  yield,  but  the  new  best  fit  line  has 
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SCALED  SPALL  THICKNESS  (mkt 


SCALED  SNOT  DEPTH  OF  BURIAL  (■it'V3t 


Figure  11.  S (maximum  height  attained  by  spall)  observed  near 
max 

the  source  vs.  scaled  shot  depth  of  burial. 


( m)  itias 


Figure  13.  S observed  near  the  source  vs.  yield, 
max 
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a smaller  y-intercept  value.  In  Figure  6 the  spall  thickness,  corrected 
to  a scaled  depth  of  burial  of  HALFBEAK  (spall  thickness  X scaled  depth  of 
burial  /122m  kt  is  plotted  versus  yield.  The  adjusted  data  show  the  same 

scatter  as  in  Figure  5.  Figure  7 shows  scaled  spall  thickness  (thickness/ 
yield^^)  versus  scaled  shot  depth  of  burial  (depth/yield^^)  . Viecelli 
found  an  inverse  relationship  between  scaled  spall  thickness  and  scaled  depth 
of  burial,  but  the  addition  of  new  data  points  shows  a much  poorer  fit  to  the 

best  fit  line  through  the  entire  data  set.  A plot  of  spall  thickness/shot  ] 

depth  versus  scaled  shot  depth  of  burial  showed  no  relationship  between  those 

two  variables.  In  Figure  8 spall  thickness  is  plotted  versus  shot  depth  of 

burial.  The  line  through  the  data  has  a slope  of  . 10±.03  and  spall  thickness 

is  observed  to  range  from  5 to  20  percent  of  shot  depth. 

In  Figure  9 the  radius  of  spallation  is  plotted  versus  yield.  Viecelli' s 
best  fit  line  is  a lower  bound  to  the  new  set  of  data.  A line  of  similar 
slope  has  been  drawn  through  the  new  set  of  data.  Both  lines  have  a slope 
of  1/3,  suggesting  the  same  dependence  on  yield  as  the  thickness  of  the 
spall.  Figure  10  shows  the  scaled  spall  radius  (radius/yield"*^)  versus 
the  scaled  shot  depth  of  burial.  No  line  could  be  found  to  fit  the  complete 
data  sot.  Viecelli' s seven  explosions  all  had  scaled  shot  depth  of  burial 
values  in  the  limited  range  95  to  105,  but  Viecelli’s  line  is  merely  a lower 
bound  for  our  data  set,  which  contains  a greater  range  of  scaled  shot  depth. 

Great  variation  exists  in  scaled  spall  radius,  even  within  a particular 
region. 

S versus  scaled  shot  depth  of  burial  is  shown  in  Figures  11  and  12. 
max 

In  Figure  11  S was  observed  near  the  source  and  in  Figure  12  S was  i 

max  ^ max 

calculated  from  observed  maximum  At.  Viecelli 's  data  exhibits  an  inverse 

I 

fourth  power  dependence  on  depth,  but  the  complete  set  of  data  in  this  study  ] 

shows  an  inverse  second  power  dependence  on  depth  and  much  more  scatter  than 

Viecelli's  data  set.  In  Figure  13  S observed  near  the  source  is  plotted  ' 

max  1 

versus  yield.  The  best-fit  line  has  a slope  of  .40±.09. 

Combining  the  results  of  Figures  5 and  9,  Figure  14  shows  spall  radius 
versus  spall  thickness.  The  data  show  a direct  relationship,  but  the  slope 
of  a line  through  the  data  is  poorly  defined.  Clearly,  data  in  Figures  5 
and  9 show  that  both  spall  depth  and  radius  are  dependent  on  yield. 


Lateral  and  vertical  changes  in  lithology  probably  are  the  major  factors  that 
affect  the  relationship  between  spall  radius  and  spall  thickness. 


Possible  relationships  between  yield  (Y),  shot  depth  (h) , spall  thick- 
ness (D) , radial  extent  of  spall  (rQ) , and  maximum  height  the  spall  is  thrown 

(S  ) were  investigated  using  a multivariate  general  linear  hypothesis 
max 

program.  The  following  relationships  were  found  between  the  variables: 

log (D)  = 1.35  + 0. 351og (Y)  - 0.071og(h) 

log(rQ)  = 2.84  + O.371og  (Y)  - 0.161og(h) 

log (S  ) = 0.57  + 0. 521og  (Y)  - 0.501og(h) 
max 

The  effect  of  shot  depth  is  significantly  different  from  zero  only  at  the  50 
percent  level,  implying  that  shot  depth  is  not  a significant  variable.  Shot 
depth  probably  is  a significant  variable  but  this  cannot  be  shown  because  of 
the  strong  correlation  of  depth  with  yield  and  the  small  number  of  observations 
(12  cases).  Both  log  (D)  and  log  (rQ)  are  approximately  proportional  to  (1/3) 
log  (Y) , as  seen  in  Figures  5 and  9. 
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ESTIMATES  OF  TELESEISMIC  pP  AND  pspALL  ARRIVAL  TIMES 
FROM  NEAR-FIELD  DATA 


Springer  (1974)  estimated  Pcn.rT-P  times  for  explosions  at  teleseismic 

o PALL 

distances  using  near-field  surface  accelerometer  data  that  show  t . the 

up 

uphole  travel-time  above  the  source,  and  At,  the  difference  between  the  spall 
arrival  time  and  the  direct  ray.  Springer  assumed  that  the  spall  signal  on 
the  surface  accelerometer  data  is  due  to  the  first,  and  deepest,  spall 
closure.  At  is  a function  of  distance  from  the  source  and  it  is  generally 
largest  over  the  source  where  the  free-fall  time  of  spall  is  largest . (Figures 
2 and  1) . pP-P  at  teleseismic  distances  can  be  approximated  by  the  formula 

pP-P  = 2t  cos0 
up 

where  e is  the  take-off  angle  at  the  source.  P___  ,-P  can  be  approximated 

SPALL 

by  the  formula 

= At  t 2t  - . 2t 
SPALL  — up  up 


(Springer  1974).  Where  no  information  on  spall  depth  is  available,  the 
travel-time  between  the  spall  depth  and  the  surface  is  assumed  to  be  20  per- 
cent of  the  uphole  time,  which  accounts  for  the  term  . 2t 

up 

Table  II  lists  pP-P  and  P<.PAT T ~P  times  for  47  explosions.  pP-P  times 
are  from  already  published  values  of  t . An  epicentral  distance  of  50 
degrees  was  used  in  the  calculations.  r -P  times  are  from  published 

near- f ield  surface  accelerometer  measurements,  particle  velocity  tables  of 
near-field  data  provided  by  Sandia  Laboratories,  and  Springer's  (1974)  study. 
The  At  values  that  Springer  used  were  near  the  shot  hole,  where  energy  was 
assumed  to  be  at  a maximum.  Since  this  is  also  the  maximum  value  of  At 
found.  Springer's  PCD.rT~P  times  aie  maximum  values.  As  demonstrated  in 

brALL 

the  next  section,  the  velocity  of  the  slab  at  impact  and,  therefore,  the 
relative  energy  of  spallation  varies  with  distance  from  the  shot  hole.  The 
impact  of  the  disc  occurs  over  a time  range  of  a few  tenths  of  a second  for 
most  explosions  and  the  maximum  energy  as  a function  of  radius  occurs  near 

Springer,  D. , 1974.  Secondary  sources  of  seismic  waves  from  underground 
nuclear  explosions,  Bull.  Seism.  Soc.  Amer.,  64 . 581. 
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TABLE  II 
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Telesexsmic  pP-P  and  P -P  Estimates. 

SPALL 


Shot 

Date 

Location 
o o 

°N  W 

Yield 

<kt) 

Depth 

of 

Burial 

(m) 

pP-P (sec) 

Acceler- 

ometers 

Particle 

Velocit les 

Springer 

RAINIER 

09-19-57 

37 

116 

1.7 

274 

.27 

.57 

.75 

ANTLER 

09-15-61 

37 

116 

2.6 

402 

.27 

.57 

SHREW 

09-16-61 

37 

116 

0-20 

98 

.11 

. 26 

MINK 

10-29-61 

37 

116 

0-20 

192 

.22 

.46 

FISHER 

12-03-61 

37 

116 

1 3 

364 

. 36 

.75 

GNOME 

12-10-61 

32 

104 

3.1 

362 

.29 

.61 

. 56 

1.32 

RINGTAIL 

12-17-61 

37 

116 

0-20 

36  3 

. 36 

.67 

STOAT 

01-09-62 

37 

116 

4.7 

302 

. 33 

.64 

DORMOUSE 

01-30-62 

37 

116 

0-20 

363 

. 36 

1.13 

ARMADILLO 

02-09-62 

37 

116 

6.5 

240 

.29 

1 .20 

HARD HAT 

02-15-6 2 

37 

116 

5.9 

287 

.12 

.44 

CHINCHILLA  I 

02-19-6 2 

37 

116 

1.8 

150 

.17 

.51 

PLATYPUS 

02-24-62 

37 

116 

0-20 

58 

.07 

.15 

PAMPAS 

03-01-62 

37 

116 

0-20 

363 

. 36 

.79 

ERMINE 

03-06-62 

37 

116 

0-20 

73 

.08 

.15 

HOGNOSE 

03-15-62 

37 

116 

0-20 

239 

.28 

.55 

CHINCILLA  II 

03-31-62 

37 

116 

0-20 

137 

.16 

.81 

DORMOUSE 

PRIME 

04-05-62 

37 

116 

n 

261 

.29 

.76 

AARDVARK 

05-12-62 

37 

116 

36 

4 34 

.27 

.80 

RACCOON 

06-01-62 

37 

116 

0-20 

164 

.19 

.47 

PACKRAT 

06-06-62 

37 

116 

0-20 

262 

.29 

1 .02 

DAMAN  I 

06-21-62 

37 

116 

0-2U 

260 

.29 

.72 

HAYMAKER 

06-27-62 

37 

116 

67 

408 

.40 

.86 

CLEARWATER 

10-16-63 

37 

116 

200-1000 

544 

.46 

.96 

SHOAL 

10-26-63 

39 

118 

12.5 

367 

.17 

.61 

.56 

.90 

SALMON 

10-22-64 

31 

90 

5.3 

828 

.47 

.86 

.86 

HANDCAR 

11-05-64 

37 

116 

12 

402 

. 30 

.85 

.92 

MUDPACK 

12-16-64 

37 

116 

2.7 

152 

.10 

.83 

MERLIN 

02-16-65 

37 

116 

10 

296 

.11 

.49 

.52 

DILUTED 

WATERS 

06-16-65 

37 

116 

0-20 

195 

.23 

1.06 

LONGSHOT 

10-29-65 

51 

179E 

85 

701 

. 37 

1 . 34 

1.76 

PINSTRIPE 

04-25-66 

37 

116 

0-20 

296 

.35 

1 .04 

DISCUS 

THROWER 

05-27-66 

37 

116 

21 

337 

.33 

2.06 

PILEDRIVER 

06-02-66 

37 

116 

%56 

463 

.19 

1.71 

1.32 

HALFBEAK 

06-30-66 

37 

116 

300 

820 

.58 

1.44 

2.93 

NEW  POINT 

12-13-66 

37 

116 

0-20 

244 

.28 

1.05 

AGILE 

02-23-67 

37 

116 

20-200 

732 

.67 

1.38 

COMMODORE 

05-20-67 

37 

116 

250 

746 

.68 

2.50 

SCOTCH 

05-23-67 

37 

116 

150 

977 

.75 

1.44 

1.73 

DOORMIST 

08-31-67 

37 

116 

0-20 

446 

.27 

.66 

LANPHER 

10-18-67 

37 

116 

20-200 

714 

.71 

2.17 

GASBUGGY 

12-10-67 

37 

107 

29 

1292 

.71 

1 .02 

1 .02 

1.17 

FAULTLESS 

01-19-68 

39 

116 

200-1000 

975 

.66 

1.26 

2.17 

KNOX 

02-21-68 

37 

116 

20-200 

645 

.55 

1.30 

2.47 

BOXCAR 

04-26-68 

37 

116 

1200 

1158 

.79 

3.18 

3.07 

CYPRESS 

02-12-69 

37 

116 

0-20 

411 

.27 

.60 

MILROW 

10-02-69 

51 

179E 

5.1000 

1219 

.61 

1.26 

1.30 

2.49 

r i 

■ 

the  shot  hole,  but  generally  not  at  the  accelerometer  closest  to  the  shot 
hole.  For  each  explosion,  the  impact  energy  was  calculated  at  each  acceler- 
ometer site  for  an  annulus  of  1 cm  width  and  a thickness  equal  to  the  spall  , 

thickness.  Where  the  energy  was  highest  the  value  of  At  at  that  acceler- 
ometer was  used  in  the  calculation  of  teleseismic  Pcn.Tr”P>  Because  the 

b r All  la 

source  of  the  Pgp^L  s i dual  extends  over  a fraction  of  a second  or  more, 

most  workers  thought  that  identification  of  the  signal  depended  upon  the 

arrival  time  of  the  maximum  energy  and  not  the  arrival  time  of  the  earliest 

energy,  which  may  be  too  small  to  be  observed  at  teleseismic  distances. 

Note  that  over  half  of  the  P_...I-P  times  are  approximately  one  second  or 

SPALL 

less.  The  separation  times  between  P,  pP  and  P may  be  too  small  1 

br  ALL 

^ these  cases  to  be  resolved  by  conventional  techniques.  Since  the  deeper 

events  would  have  greater  times  of  separation  between  P,  pP,  and  P 

b I Alill  f 

the  best  chances  for  detection  of  pP  and  PgpAI j are  for  the  deeper  events 
(depth  of  burial  greater  than  500  m at  NTS) . 
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SPALL  ENERGY  AND  ITS  EFFECTS  ON  m AND  M 

b s 

Energy  (ergs)  as  a function  of  range  was  calculated  for  all  explo- 
sions using  the  formula 

ENERGY (r)  = 27frDAr(52gh  (r) 

where  r is  the  horizontal  range,  D is  the  spall  thickness,  6 is  the  over- 
burden density,  g is  gravitational  acceleration,  and  h(r)  is  the  height  to 
which  the  spall  is  thrown.  For  each  explosion  the  impact  of  an  annulus  of 
one  cm  diameter  ( Ar  equals  one  cm)  was  calculated  at  each  accelerometer  site. 
Figure  15  shows  estimated  energy  as  a function  of  range  for  five  explosions. 
These  five  explosions,  shown  in  a later  section,  have  relatively  hiqh 
spall  energy  compared  to  the  energy  of  the  explosion.  No  relationship 
exists  between  the  range  of  maximum  energy  and  spall  thickness,  spall 
radius,  or  geographic  location.  However,  there  was  direct  relationship 
between  yield  and  tne  range  of  maximum  spall  energy  (Figure  16). 

For  the  entire  set  of  47  explosions,  the  radius  of  maximum  spall  energy 
ranged  from  9 to  100  percent  of  the  maximum  radius  of  spall. 

Figure  17  shows  the  energy  values  shown  in  Figure  15  plotted  as  a 
function  of  time  after  the  origin  time  of  the  explosion.  Each  vertical  line 
represents  the  energy  of  impact  at  the  radius  of  an  accelerometer  site  for  an 
annulus  of  one  cm  diameter.  The  data  has  been  sampled  discontinuously  at  each 
of  the  near-field  accelerometer  sites  and  the  dashed  lines  between  energy 
measurements  suggest  continuous  production  of  energy  as  the  slab  impacts 
over  a time  interval  of  less  than  a second  for  these  five  explosions.  The 
teleseismic  spall  signal  would  be  emergent  with  smaller  energies  from 
the  closures  at  large  distances  from  the  shot  hole  that  arrive  befcre  the 
higher  energies  from  closures  at  closer  distances  to  the  shot  hole  where 
• free-fall  time  is  greatest. 

The  author  wanted  to  examine  the  relative  energy  of  spallation  compared 
to  the  explosion's  energy.  From  a knowledge  of  the  mass  of  the  falling  slab 
and  the  height  to  which  it  is  thrown,  the  total  energy  in  a falling  slab  can 
be  calculated.  Tables  that  sandia  Laboratories  provided  gave  the  height  the 
slab  was  thrown  as  h^  atdistance  r^.  The  slab  is  divided  into  annuli 
centered  about  the  explosion.  The  radii  from  the  explosion  to  the  outer 
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and  inner  boundaries  of  the  annuli  would  be  the  radii  where  the  height 
of  the  thrown  slab  was  measured.  By  summing  the  individual  contributions 
from  each  falling  annulus  the  total  energy  would  be  found 


2 

ENERGY  (ergs)  = Y.  2trr  iDAri<SVQ 

i=l 

where  N = number  of  annuli  or  number  of  observation  sites 

l r.  = the  distance  to  the  center  of  each  annulus, 

| 1 

r . = R . + K . 

l — l l-l 

D = thickness  of  annuli  (spall  thickness) 

Ar  = width  of  annuli,  r.  = P.  - R.  , 
t li  l-l 

■S  - density  of  annuli 

V = velocity  of  slab  at  impact, 

V“  = 2gh(r)  and  h(r)  = H.  + H. 

0 — l l + l 

2 

If  the  spall  thickness  was  not  available  from  vertical  accelerometer  arrays, 
or  observed  parting  in  the  boreholes,  then  spall  thickness  was  assumed  to  be 
IQ  percent  of  the  shot  burial  depth,  a figure  suggested  by  the  relationship 
between  spall  thickness  and  shot  depth  of  burial  in  Figure  8.  Table 
III  shows  the  spall  and  explosion  energies  for  the  explosions  with  known 
yields.  Explosion  energies  were  found  in  Ferret  (1872)  or  were  calculated  | 

with  the  efficiencies  for  various  mediums  found  in  Ferret's  study  according 
to  the  formula 

19 

ENERGY(ergs)  = Yield(kt)  x Efficiency  x 4.18x10  (1) 

To  examine  the  amount  of  energy  that  spall  generated  compared  to  the  energy 
of  the  explosion,  the  ratio  of  these  two  values  for  each  explosion  was  plotted 
versus  yield  (Figure  18).  The  explosion  had  spall  energies  that  were  from  1 to 
J22  percent  of  the  explosion  energies.  Obviously,  errors  exist  in  the  assump- 
tions used  to  calculate  the  spall  energy  or  the  radiated  energy  since  the  spall 
energy  cannot  exceed  that  fraction  of  the  explosion  energy  traveling  upward  into 

Ferret,  W.,  1972.  Seismic-source  energies  of  underground  nuclear  explosions. 

Bull.  Seism.  Soc . Amor.,  62 . 763. 
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TABU-:  III. 


Spall  Energies 


Explosion  Energy  (10*°orgs)  E, 


Explosion 

Spa  1 1 
Energy 
(10  ergs) 

Ferret 

Ca leu lat od* 

E 

EXPLOSION 

M 

s 

(von  Seggern) 

RAINIER 

.027 

1.5 

. 18 

ANTLER 

.012 

2 . 2 

. 06 

FISHER 

.002 

.68 

.03 

GNOME 

.214 

12. 

. 18 

ARMADILLO 

.0135 

27.2 

.01 

HARDHAT 

.0245 

8.4 

.03 

3.01 

CHINCHILLA  I 

.0050 

7.5 

.01 

AARPVARE 

. 666 

30.1 

. 22 

3.56 

HAYMAKER 

. 324 

2.7 

1 .20 

3.02 

SHOAL 

.841 

13. 

.65 

SALMON 

.07  5 

7.6 

. 10 

HANDCAR 

.444 

11. 

.40 

3.03 

Mi  DDACK 

.008 

. 12 

.66 

1.76 

merlin 

.040 

. 17 

2.36 

I .ON  '.SHOT 

7.77 

71. 

1.09 

DISCI'S  THROWER 

.8  36 

2.6 

3.22 

3.05 

FTI.EPKIVEK 

3.79 

47. 

.81 

4.05 

HALFBEAK 

41.1 

251  . 

1.64 

5.06 

COMMODORE 

7.01 

209. 

.34 

4.59 

SCOTCH 

8.00 

125. 

.64 

4.51 

GASBUGGY 

.850 

24. 

.35 

BOXCAR 

226. 

1400. 

1 .61 

5.56 

Mil  ROW 

54.9 

8 36. 

.72 

5.0 

19 

* Energy  - Yield  x Efficiency  x 4.18  x 10 


Figure  18.  Spall  energy  /explosion  energy  vs.  yield. 


the  spall  region.  Assuming  a spherical  explosion  source,  the  fraction  of 
the  explosion  energy  reaching  the  spall  region  depends  upon  the  depth  of 
burial  and  the  spall  radius.  For  explosions  studied  here  that  fraction 
ranged  from  1 to  37  percent. 

Spall  energy  was  calculated  for  42  NTS  explosions  located  below  the 

water  table  that  lacked  near-field  data.  Thickness  of  the  spall  layer  was 

estimated  from  the  relationship  between  yield  and  spall  thickness 

found  in  Figure  5.  Similarly,  the  maximum  radius  of  spall  was  found  from 

Figure  9 and  the  maximum  height  attained  by  spall  was  found  from  Figure  11. 

The  height  spall  was  thrown  was  assumed  to  vary  inversely  as  the  square  of 

the  horizontal  distance  from  the  shot  hole,  which  approximates  the  conical 

shape  of  soall  height  observed  by  other  investigators  (Figure  31).  Figure 

l'1  shows  the  M residual  (M  - log  yield  + constant.)  for  these  explosions 

(von  Seggern  and  Gurski,  197(3)  versus  spall  enerqy/oxplosion  energy. 

No  relationship  exists  between  M residual  and  spall  energy/explosion  energy, 

s 

a situation  suggesting  that  spall  has  no  effect  on  M . In  addition,  low 
strain  release  events  (Masse,  unpublished  data)  do  not  show  a clearer 
relationship  titan  do  average  strain  release  events,  suggesting  that  strain 
release  is  not  hiding  the  effect  of  spall  energy. 

The  efficiency  that.  Ferret  (1972)  calculated  is  a source  of  error 
in  the  energy  calculations.  Table  XV  shows  the  explosion  energies  that 
Ferret  derived  from  near-field  data  and  the  explosion  energies  calculated 
for  the  same  explosions  using  formula  (1).  The  range  of  variation  between 
the  two  energy  values  for  each  explosion,  and  hence  the  error  in  estimating 
explosion  energies  using  Perret's  efficiencies,  ranged  from  a factor  of  1.0 
to  2.7. 

Another  source  of  error  in  the  energy  calculations  is  the  thickness  of 
the  spall  layer.  In  a few  cases  partings  were  observed  in  boreholes.  However, 
there  were  several  spall  partings  and  the  one  responsible  for  the  spall  signal 
on  the  near-field  surface  accelerometer  data  must  be  determined.  Because  most 

von  Seggern,  D.,  and  J.  Gurski,  1976.  M versus  yield  of  underground  nuclear 
explosions  at  the  Nevada  Test  Sites,  SDAC-TR-76-11 , Teledyne  Geotech, 
Alexandria,  Virginia  22314. 
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Figure  19.  residual  vs.  spall  energy/explosion  energy  for  42  NTS 

explosions  located  below  the  water  table. 
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TABLE  IV 


Explosion  energies  from  Perret  (1972)  and  from 
Ferret's  efficiency  estimates  for  the  explosion  medium. 


Shot 

Enerqy 

Perret 

(10  ergs) 

Yield  x efficiency 
, x 41.8 

Dif fer- 
erence 

..  Lil 

Yield 

Pff  i r-i  enny 

RAINIER 

1.5 

1.4  - 2.1 

0-40 

1.7 

2-3 

FISHER 

.68 

.27  - .82 

0-60 

13 

.05  - .15 

GNOME 

12. 

3.9 

68 

3.1 

3 

HARDHAT 

8.4 

4.9 

42 

5.9 

2 

HAYMAKER 

2.7 

1.40  - 4.2 

0-56 

67 

.05  - .15 

SHOAL 

13. 

10.4 

20 

12.5 

2 

SALMON 

7.6 

6.6 

13 

5.3 

3 

HANDCAR 

11. 

10.0 

9 

12 

2 

MUDPACK 

. .12 

.11  - .34 

0-183 

2.7 

.10  - .30 

MERLIN 

.17 

.21  - .63 

24-271 

10 

.05  - .15 

DISCUS  THROWER  2.6 

.89  - 2.6 

0-66 

21 

.10  • .30 

PILEDRIVER 

47. 

47. 

0 

56 

2 

GAS BUGGY 

24. 

24. 

0 

29 

2 

BOXCAR 

1400. 

1003. 

28 

1200 

2 
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subsurface  displacements  are  attributed  to  parting  at  one  depth,  the  energy 
for  the  falling  disc  varies  directly  with  the  thickness  of  the  falling  disc. 


In  some  cases  spall  depth  was  determined  from  observations  of  the 

arrival  times  of  the  phase  t^  and  the  phase  reflected  off  the  spall  layer  t^. 

Figure  2 illustrates  the  ray  paths  of  these  two  phases.  This  method  of  \ 

determining  spall  depth  is  suspect  because  identification  of  the  presumed 

t-  arrival  is  uncertain.  If  the  sources  of  the  spall  signals  were  close 
3 1 
to  the  surface,  then  the  t^  arrival  could  not  be  separated  from  the  t^  arrival, 

and  the  total  amount  of  energy  in  such  a thin  falling  disc  would  be  j 

small  compared  to  the  explosion  energy.  i 

I 

1 

In  cases  where  subsurface  vertical  accelerometer  arrays  were  used  j 

I 

to  determine  spall  depth,  the  spall  depth  was  located  between  two  sensor  , 

depths  with  the  spoil  signal  showing  moveout  in  time  above  the  presumed  spall 
depth.  Figures  20  through  23  show  examples  of  near-field  subsurface  vertical 
accelerometer  array  data.  In  each  case  note  the  arrival  of  the  initial  accel-  . 

eration  pulse,  travelling  upward  at  the  compressional  velocity  of  the  medium. 

Dashed  lines  indicate  expected  return  of  the  reflected  phase  pP.  In  all 
cases  there  is  little  pP  energy  leaving  the  source  region  compared  to  the 
initial  signal,  suggesting  that  for  these  explosions  pP  is  of  small  amplitude 
compared  to  P at  teleseismic  distances.  Note  that  pP  has  travelled  farther 
than  P and  is,  therefore,  expected  to  be  smaller  because  of  spherical 

divergence.  The  next  larqe  arrival  following  the  initial  pulse  is  the 
spall  signal.  GNOME  and  SLAMON  show  several  spall  signals.  For  RAINIER, 
assuming  that  all  spall  occurs  above  these  depths,  the  data  (Figure  20) 
shows  no  large  spall  pulses  at  137-,  160-  and  162-  meter  gauges.  Similarly, 
spall  is  assumed  to  occur  above  91  meters  for  SALMON  (Figure  22)  and  above 
1000  feet  for  MILROW  (Figure  23).  Spall  closure  signals  in  Figures  20 
through  23  propagate  upward  at  a speed  slower  than  the  compression  velocity 
of  the  medium.  Thus,  spalling  is  presumed  to  occur  over  a range  of  depths. 

Deeper  formations  close  before  shallower  layers  and  accelerometer  data  show 

each  spall  phase  occurs  in  the  first  cycle  caused  by  the  spall  closure  closest 

to  that  gauge.  The  total  energy  in  spall  closure  probably  results  from 

summing  at  least  several  spall  closures  and  not  from  initial  spall  impact.  < 
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Figure  20.  Subsurface  vertical  accelerometer  array  data 

for  RAINIER  (Eisler,  Chilton,  and  Sauer,  1966) . 
Arrows  indicate  the  direction  of  flow  of  seismic 
energy. 
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Figure  21.  Subsurface  vertical  accelerometer  array  data 
for  GNOME  (Weart,  1962) . Lower  gauges  suffer 
from  early  cable  failure.  Arrows  indicate  the 
direction  of  flow  of  seismic  energy. 
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Figure  23.  Subsurface  vertical  accelerometer  array  data 

for  MILROW  (Perret  and  Breding,  1972) . Horizontal 
range  = 299  feet. 
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If  spall  energy  is  large  enough  to  be  observed  at  teleseismic  distances,  then 
after  closing  each  layer  the  spall  signal  should  be  propagating  downward  at 
the  compressional  velocity  of  the  medium.  Dashed  lines  in  Figures  20  through 
23  indicate  expected  downward  flow  of  spall  energy  but  none  of  these  expected 
signals  are  observed.  This  situation  exists  because  the  spall  signal  for 
each  layer  closure  is  extended  over  a few  tenths  of  a second,  and  each 
boundary  is  closing  first  at  large  ranges  from  the  source  and  then  near  the 
shot  hole  where  free-fall  times  are  largest.  If  spall  energy  is  significant, 
then  spall  signals  observed  at  the  surface  should  show  extended  periods  of 
high  amplitudes,  but  instead  each  signal  is  largely  one  pulse.  Sources  of 
spall  signals  are  probably  relatively  low  in  energy  and  located  close  to  the 
gauges.  The  spall  signals  only  seem  significant  on  surface  records  when 
compared  to  the  initial  acceleration  pulse  because  the  initial  pulse  has 
travelled  a longer  distance  to  the  guage  and  is  more  highly  attenuated  than 
the  spall  signal.  In  fact,  the  subsurface  vertical  accelerometer  array  data 
show  little  if  any  spall  energy  leaving  the  source  region.  Thus  spall 
appears  to  have  little  effect  on  and  at  teleseismic  distances. 

Subsurface  observations  of  spall  signals  should  be  interpreted  as 
observations  of  spall  closures  closest  to  that  gauge  and  not  observations 
of  the  first,  and  deepest,  spall  closures.  Thus,  the  depths  and  vertical 
displacements  of  spall  openings  can  only  be  determined  with  vertical  sub- 
surface arrays  of  accelerometers.  Consequently,  calculations  of  total 
energy  in  the  closure  of  all  spall  layers  depend  upon  subsurface  data. 

P -P  times  based  on  surface  data  are  valid  only  for  closure  of  the 

jPALjIj 

shallowest  layer,  which  may  not  be  the  spall  closure  containing  the  most 
impact  energy.  The  total  height  the  ground  surface  is  thrown  is,  then,  the  sum 
of  multiple  spall  openings  and  it  is  not  due  to  the  first,  and  deepest, 
spall  opening. 


CONCLUSIONS 


1 


Sandia  Laboratories  near-field  observations  of  spall  on  particle  velocity 
data  were  used  to  estimate  possible  relationships  between  radius  of 
spallation,  thickness  of  spall,  the  maximum  height  attained  by  spall,  shot 
depth  of  burial,  and  yield  for  47  explosions.  The  following  relationships 
were  found 

log  (spall  thickness)  a 1/3  log  (yield) 
log  (spall  radius)  a( 1/3) log  (yield) 
log  (maximum  spall  height)  a .40  log  (yield) 
spall  thickness  a -10  shot  depth  of  burial 

Maximum  spall  height  is  inversely  proportional  to  the  scaled  shot  depth  of 

burial  but  a best  fit  line  throuqh  the  data  is  poorly  defined.  Most  scatter 

of  data  points  in  all  comparisons  is  due  probably  to  varying  lithologies  ir 

in  the  source  regions.  Estimates  of  pP-p  and  P , - P times  at  teleseis- 

SpaLL 

mic  distances  based  upon  surface  accelerometer  data  are  given  in  Table  II. 

The  spall  arrival  times  are  from  0.2  to  3.3  seconds  after  the  main  phase. 


"sing  near-field  displacement  data,  energies  of  the  falling  disc  were 
calculated  for  each  explosion  and  then  compared  to  explosion  energy  that 
Perret  (1972)  calculated,  or  derived  from  the  efficiencies  found  by  Perret 
for  various  mediums.  The  explosions  had  spall  energies  ranging  from  1 to 
322  percent  of  the  explosion  energies.  Thus,  errors  exist  in  assumptions 
used  to  calculate  spall  energy  and/or  explosion  energy  because  spall  energy 
cannot  exceed  that  fraction  of  the  explosion  energy  traveling  upward  into  the 
spall  region;  for  explosions  in  this  study  that  fraction  ranged  from  1 to  37 
percent.  Relative  spall  energies  compared  to  explosion  energies  for  42  NTS 
explosions  located  below  the  water  table  showed  no  relationship  to  M^  or 
strain  release. 


4 


Subsurface  vertical  accelerometer  data  suggests  that  spall  occurs  over 
a range  of  depths,  the  deeper  formations  closing  before  the  shallower  layers. 
The  large  observed  pulse  of  each  spall  closing  is  caused  by  the  spall  closure 
closest  to  that  gauge.  The  surface  accelerometer  records  indicate  large 
spali  energy,  but  the  initial  pulse  has  traveled  a longer  distance  to  the 
gauge  and  it  has  been  more  highly  attenuated  than  the  spall  signal.  The 
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principal  error  in  the  energy  calculations  was  the  assumption  that  one  spall 
parting  was  largely  responsible  for  observed  surface  displacements.  Total 
energy,  the  sum  of  many  spall  closures,  is  probably  much  smaller  than  the 
spall  energies  calculated  assuming  that  the  first,  and  deepest,  spall 
closure  contained  the  most  energy.  Surface  observations  of  spall  have 
often  been  interpreted  incorrectly  by  other  investigators.  Spall  signals 
observed  at  the  surface  are  observations  of  the  spall  closure  closest  to  that 
gauge,  not  necessarily  observations  of  the  first,  and  deepest,  spall  closure. 
P — P times  calculated  using  surface  data  are  for  spall  closest  to  the 

b r 

surface,  which  is  the  last  spiall  to  close  and  may  not  be  the  spall  closure 
containing  the  most  energy.  Determining  depths  and  vertical  displacements 
of  spall  openings,  and  hence  calculation  of  total  energy  in  spall,  depends 
upon  observations  from  vertical  subsurface  accelerometer  arrays.  The 
total  height  the  ground  surface  is  thrown  is  the  sum  of  multiple  spall 
openings  and  it  is  not  due  to  the  first,  and  deepest,  spall  opening.  The 
subsurface  accelerometer  data  show  little,  if  any,  spall  energy  leaving 
the  source  region,  suggesting  that  spall  has  little  effect  on  m^  and  at 
teleseismic  distances.  Further  investigations  should  accumulate  more  sub- 
surface accelerometer  records  and  filter  the  data  to  enhance  the  lower 
frequencies  which  are  more  important  in  teleseismic  analysis  than  the  high 
frequency  results  shown  in  the  literature. 
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Surface  Particle  Velocity  Data 
(Sandia  Laboratories) 
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Explosion 


Surface  Particle  Velocity  Data 
(Sandia  Laboratories) 

Horizontal  . 5At  Peak  Displacement 

Range  (m)  (sec)  (cm) 


RAINIER 

30 

.25 

10.62 

458 

.14 

1.60 

ANTLER 

149 

.20 

3.56 

301 

.14 

9.40 

SHREW 

0 

.07 

.71 

MINK 

9 

.11 

2.54 

FISHER 

23 

.18 

9.91 

GNOME 

402 

.12 

13.72 

805 

.11 

3.05 

RINGTAIL 

23 

.14 

3.81 

STOAT 

61 

.14 

4.32 

DORMOUSE 

46 

.20 

11.68 

198 

.37 

14.22 

366 

.32 

7.62 

ARMADILLO 

46 

.22 

21.34 

122 

.44 

8.64 

244 

.15 

3.56 

HARDHAT 

152 

.31 

5.59 

305 

.15 

15.49 

CHINCHILLA  I 

15 

.40 

33.02 

84 

.16 

7.37 

168 

.16 

1.78 

335 

.12 

.51 

PLATYPUS 

15 

.05 

.86 

32 

.04 

.51 

PAMPAS 

15 

.20 

17.02 

ERMINE 

15 

.03 

.43 

A - 1 
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Surface  Particle  Velocity  Data 
(Sandia  Laboratories) 


Explosion 

Horizontal 

Ranqe  (m) 

. 5At 
( sec) 

Peak  Displacement 
(cm) 

HOGNOSE 

9 

.17 

8.89 

183 

.45 

5.33 

274 

. 32 

8.89 

CHINCHILLA  II 

15 

.20 

16.00 

69 

. 32 

7.11 

137 

.09 

2.54 

DORMOUSE  PRIME 

15 

.22 

20.32 

AARDVARK 

15 

.45 

55.88 

220 

.26 

33.02 

440 

.25 

19.56 

880 

.18 

8.13 

RACCOON 

15 

.13 

12.95 

PACKRAT 

15 

.35 

12.70 

DAMAN  I 

15 

.42 

25.40 

132 

.20 

16.51 

HAYMAKER 

15 

.59 

43.18 

198 

.25 

12.95 

412 

.26 

11.18 

824 

.21 

6.35 

CLEARWATER 

0 

.29 

60.96 

211 

. 34 

30.48 

364 

.28 

25.15 

668 

.23 

16.51 

SHOAL 

92 

.17 

57.40 

243 

.15 

21.34 

396 

.19 

32.00 

596 

.13 

4.57 

SALMON 

24 

.17 

34.04 

149 

.20 

25.91 

168 

.17 

25.91 

320 

.11 

9.14 
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Explosion 


Surface  Particle  Velocity  Data 
(Sandia  Laboratories) 

Horizontal  . 5At 

Range  (m)  (sec) 


Peak  Displacement 
(cm) 


HANDCAR 


MUDPACK 


MERLIN 


DILUTED  WATERS 


LONGSHOT 


PINSTRIPE 


DISCUS  THROWER 


15 

.31 

34.80 

122 

.32 

27.43 

213 

.26 

20.32 

427 

.11 

6.35 

15 

.39 

39.37 

46 

.43 

25.91 

91 

.36 

13.21 

168 

.39 

4.32 

259 

.09 

1.17 

15 

.25 

11.61 

30 

.23 

13.21 

46 

.22 

11.91 

76 

.22 

12.70 

107 

.22 

11.02 

152 

.20 

10.92 

305 

.14 

4.93 

457 

.11 

1.70 

15 

.61 

73.66 

76 

.40 

78.74 

152 

.78 

38.10 

229 

.61 

16.26 

500 

.50 

7.37 

6 

.72 

224.28 

183 

.56 

90.42 

335 

.47 

91.19 

792 

.20 

30.99 

1219 

.20 

12.70 

1798 

.12 

5.08 

0 

.55 

61.47 

91 

.48 

61.21 

152 

.33 

42.93 

198 

.11 

17.27 

15 

.85 

213.36 

46 

.78 

172.72 

122 

.85 

132.08 

488 

.58 

4.57 

1341 

.21 

1.27 
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APPENDIX  Continued 

Surface  Particle  Velocity  Data 
(Sandia  Laboratories) 

Explosion  Horizontal  . 5At  Peak  Displacement 

Range  (m)  (sec)  (cm) 


PILEDRIVER 

9 

.65 

27 

.73 

110 

.66 

176 

.65 

226 

.75 

306 

.60 

368 

.48 

432 

.43 

573 

.40 

744 

.21 

1372 

.13 

HALFBEAK 

15 

.91 

457 

.40 

914 

. 36 

2134 

.20 

NEW  POINT 

6 

.34 

30 

. 30 

46 

.37 

91 

.26 

198 

.14 

AGILE 

13 

.65 

91 

.64 

457 

.40 

905 

.32 

1432 

. 22 

COMMODORE 

15 

1.24 

91 

1.21 

457 

.88 

969 

.28 

SCOTCH 

22 

.44 

856 

.32 

i 206 

.31 

2291 

. 17 

DOORMIST 

0 

.35 

125 

.28 

226 

.18 

320 

.15 

233.68 

254.00  . 1 

271.78 

243.84 

248.92 

157.48 

142.24  ; 

106.68 
66.04 
22.86 

3.81 

480.06 

157.48 
35.56 
10.16 

25.65 

20.83 

36.58 

12.45 

6.60 

99.06 

77.98  i 

45.21 
36.32 
8.89 

378.46  j 

99.06 

125.48 

38.10  i 

i 

89.41 

25.91 

34.29  * 

8.13  i 

i 

73.66  i 

33.02  j 

35.56  j 

7.37  j 
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APPENDIX  Continued 


Surface  Particle  Velocity  Data 
(Sandia  Laboratories) 

Explosion  Horizontal  . 5At  Peak  Displacement 

Ranqe  (m)  (sec)  (cm) 


LANPHER 

15 

.83 

152.40 

91 

.78 

167.64 

4 57 

. 70 

63.50 

900 

. 30 

28.45 

1437 

.26 

15.00 

GASBUGGY 

30 

.20 

17.09 

130 

.20 

15.06 

461 

.17 

8.66 

802 

.12 

6.91 

FAULTLESS 

26 

.82 

259.00 

152 

.59 

108.97 

457 

.53 

124.97 

1036 

.44 

89.92 

1828 

.37 

36.07 

3647 

.27 

28.96 

KNOX 

15 

1.08 

304.80 

91 

1.02 

276.86 

457 

. 35 

58.93 

BOXCAR 

23 

.93 

327.66 

305 

1.16 

477.52 

914 

.49 

133.35 

2404 

. 36 

33.02 

3813 

.31 

15.24 

CYPRESS 

0 

.45 

34.80 

125 

.15 

24.38 

MILROW 

76 

.95 

431.80 

573 

.50 

149.86 

1225 

.37 

78.74 

1354 

. 32 

81.28 

1837 

.27 

27.94 

2010 

.22 

17.53 

